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ABSTRACT 



We report on the first simultaneous observations that cover the optical, X-ray, 
and high energy gamma-ray bands of the BL Lac object PKS 2155—304. The 
gamma-ray bands were observed for 11 days, between 25 August and 6 September 
2008 (MJD 54704-54715), jointly with the Fermi Gamma-ray Space Telescope 
and the H.E.S.S. atmospheric Cherenkov array, providing the first simultaneous 
MeV-TeV spectral energy distribution (SED) with the new generation of 7-ray 
telescopes. The ATOM telescope and the RXTE and Swift observatories provided 
optical and X-ray coverage of the low-energy component over the same time 
period. The object was close to the lowest archival X-ray and Very High Energy 
(VHE, > 100 GeV) state, whereas the optical flux was much higher. The light 
curves show relatively little (~ 30%) variability overall when compared to past 
flaring episodes, but we find a clear optical/VHE correlation and evidence for 
a correlation of the X-rays with the high energy spectral index. Contrary to 
previous observations in the flaring state, we do not find any correlation between 
the X-ray and VHE components. Although synchrotron self-Compton models are 
often invoked to explain the SEDs of BL Lac objects, the most common versions 
of these models are at odds with the correlated variability we find in the various 
bands for PKS 2155-304. 

Subject headings: Galaxies: active - BL Lacertae objects: Individual: PKS 2155— 304 
- Gamma rays: observations 



1. Introduction 

The underlying particle distributions of blazars are usually studied by matching broad- 
band observations with predictions from radiative models. Since these sources are highly 
variable, simultaneous observations are essential. The most energetic BL Lac spectra extend 
up to TeV energies, and positive detections have usually indicated flaring states. However, 
with their improved sensitivity, the new generation of Atmospheric Cherenkov Telescopes 
(ACTs), which has more than quadrupled.] the number of known extragalactic VHE sources, 



50b School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden 

1 See, e.g., the online TeVCat catalog [http: / / tevcat.uchicago.edu[ which has 22 sources at the time of the 
writing of this Letter. 
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finds a few of these sources in marginally variable states with consistent detections after short 
exposures. One of these objects, the blazar PKS 2155—304 at z = 0.116, is an ideal target for 
such studies. Crucial information is expected from the Fermi Gamma-ray Space Telescope, 
since its improved sensitivity over EGRET would constrain dramatically the existing models 
that predict a wide variety of fluxes in the 100 MeV-10 TeV energy range. Since the H.E.S.S. 
experiment detects this source in a low state within ~ 1 h, significant daily detections were 
guaranteed and the source was targeted for a 11-day multiwavelength campaign. 



2. Observations and Analysis Results 



The H.E.S.S. observations of PKS 2155-304 took place during MJD 54701-54715, for 
a total of 42.2 hours. After applying the standard H.E.S.S. data-quality selection criteria, 
an exposure of 32.9 hours live time remains (MJD 54704-54715), at a mean zenith angle 
of 18.3°. The data set has been calibrated using the standard H.E.S.S. calibration method 
(jAharonian et al.l 12004 ) . The analysis tools and the event-selection criteria used for the 
VHE analysis are presented in lAharonian et al.l (120081 ). The events have been selected us- 
ing "loose cuts" , preferred for their lower energy threshold of 200 GeV and higher 7-ray 
acceptance. A 0.2° radius circular region centered on PKS 2155—304 was defined to collect 
the on-source events. Th e background was estimated using the "Reflected Region" method 
(jAharonian et al.ll2006bl ). Those observa tions yield an ex cess of 8800 events, a signal with a 
significance of 55.7<r calculated following iLi fc Ma I (119831 ) . Using standard cuts an excess of 
3612 events wi th a significance of 68.7<r is found. An independent analysis and calibration 
( iBenbowl 120051 ) yields similar results. 



The data from the Large Area Telescope (LAT; lAtwood et al.ll2008l ) have been analyzed 
by using ScienceTools v9.7, which will be publicly available from the HEASARC in the 
future. Events having the highest probability of being photons (class 3, called "diffuse") 
and coming from zenith angles < 105° (to avoid Earth's albedo) were selected. The diffuse 
emission along the plane of the Milky Way, mainly due to cosmic-ray interactions with the 
Galactic interstellar m atter, has been modeled using the 54 59Xvarh7S model pr epared with 
the GALPROP code Jstrong. Moskalenko k Reimer I booi Istrong et all l2004h which has 
been refined with Fermi-LAT data taken during the first 3 months of operation. The extra- 
galactic diffuse emission and the residual instrumental background have been modeled as an 
isotropic power-law component and included in the fit. Photons were extracted from a region 
with 10° radius centered on the coordinates of PKS 2155—304 an d analyzed with an unbinned 
maximum likelihood technique (jCash Ill979l ; iMattox et al.lll996l ) using the Likelihood anal- 
ysis software provided by the LAT team. Because of calibration uncertainties at low energies, 
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data in the 0.2-300 GeV energy band were selected. 



A total of 75 ks of exposure was taken with RXTE, spread over 10 days coinciding 
with the H.E.S.S. observations, and a 6.4 ks exposure with Swift was ma de towards the 
end of the campaign . The data taken with the PCA (IJahoda et al.l 119951 ) and the XRT 
(IBurrows et al.l 120051 ) instruments were analyzed using the HEASOFT 6.5.1 package using 
the Guest Observer Facility recommended criteria. The XRT data were extracted from a 
56" slice, both for the source and the background. Since the rate was less than 10 Hz, no 
pile-up is expected in the Windowed Timing (WT) mode. 

During the multiwavelengt h campaign a tota l of 106 observations were taken with the 
0.8 m ATOM optical telescope (IHauser et al.ll2004l ) located on the H.E.S.S. site. Integration 
times between 60 s and 200 s in the Bessel BVR filter bands were used. Photometric accuracy 
is typically between 0.01 and 0.02 mag for BVR. 



2.1. Spectral Analyses 



The H.E.S.S. time- ayeraged photon s pectrum is derived using a forward-folding maxi- 
mum likelihood method (IPiron et al.ll200ll ). The very high energy data are well described by 
a power-law of the form dN/dE = I Q (E / E )~ r , with a differential flux at E = 350 GeV (the 
fit decorrelation energy) of Iq = 10.4±0.24 stat ± 2.08 sys x 10~ n cm -2 s _1 TeV -1 and a spectral 



index T = 3.34 ± 0.05 sta t ± 0.1 sys . As before, during non-flaring states of PKS 2155— 304, 
the spectrum, measured with limited event statistics, shows no indication of curvature. The 
spectral index is similar to that p reviously measured by H .E.S.S. when the source was at 
a comparable flux leve l, in 2003 JAharonian et al.l l2005al bh and between 2003 and 2005 
( jAharonian et al.l 120081 ). The VHE spectrum is affected by interactions with the Extragalac- 
tic Backgrou nd Light (EBL) which modifies the intrinsic shape and intensity. Using the 
P0.45 model (jAharonian et al.ll2006al ). the intrinsic spectral index is derived to be r; nt ~ 2.5. 



The average Fermi spectra over the duration of the campaign are fitted by a simple power 
law for which I„ = (2.42 ± 0.33 stat ± 0.16 sys ) x lO^cm" 2 s" 1 MeV" 1 , T = 1.81 ± 0.11 sta t ± 
0.09 sys , and Eq = 943 MeV is the energy at which the correlation between the fitted values of 
T and Iq is minimized. The total exposure is 7.7 x 10 8 cm 2 s. There is no statistical preference 
for a broken power law in this data set. The light curve derived for Fermi data between 
MJD 54682-54743 shows a similar state on average as during this campaign, so in order to 
increase the photon statistics for the spectral fits, those data were included, resulting in an 
increase of the exposure by a factor of 3.6. The longer data set is then fit by a broken power 
law spectrum, which is preferred over the single power law with a significance of 97% using 
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the likelihood ratio test. We obtain a low-energy photon index of Tl — 1.61±0.16 sta t±0.17 sys , 
a break energy of E br = 1.0 ± 0.3 GeV, a high-energy index of T H = 1.96 ± 0.08 sta t ± 0.08 sys , 
and a 0.2-300 GeV flux of (1.13 ± 0.05 stat ± 0.11 sys ) x 10" 7 cm" 2 s" 1 . The Fermi spectrum 
is consist ent with the hard photon index of 1 .71 ± 0.24 during a flaring episode detected by 
EGRET JVestrand. Stacv fc Sreekumari ll995h. but it differs from the Third EGRET Catalog 
spectrum (IHartman et al.lll999l ) where the index is 2.35 ± 0.26. 

The 4-10 keV PCA and 0.5-9 keV XRT data were analyzed simultaneously with XSPEC 
V12.4.0 JArnaudlE996l ). using a broken power law model and taking into account the un- 
certainty in the cross-calibrations, as well as the variability across the non-simultaneous 
observation s, by using a multiplicative factor fo r each instrument (fixed to 1 for the PCA 



data) as in 
of Nh = 1.48 x 10 



Falanga. Belloni fc Campanal (120061 ) . Using a fixed Galactic hydrogen column 



-20 



cm 2 , we obtain a low-energy photon index of 



2.36 ± 0.01, a 

break energy of Et, r = 4.44 ± 0.48 keV, and a high-energy index ofl^ = 2.67 ±0.01, for 
an unabsorbed 2-10 keV flux of 4.99 x 10~ n ergs cm~ 2 s" 1 , wh ich is approximately 2 times 
higher than during the 2003 campaign (lAharonian et al.ll2006al ). This is similar to the VHE 
flux increase reported above, while s till being well below the high state fluxes reported by 



Vestrand. Stacy fc Sreekumarl (119951 ). 



2.2. Light Curves 



The light curves from H.E.S.S., Fermi, RXTE and ATOM are shown in Fig. [lj where 
the H.E.S.S. runs (~ 28min) were combined to derive nightly flux values. The average 
integrated flux above 200 GeV, (5.56 ± 0.13 sta t ± l-ll sys ) x 10~ n phcm~ 2 s" 1 , corresponds 
to ~ 20% F Crab>2 ooGeV, or ~ 50% higher than the qui escent state of 2003 (lAharonian et al. 



2006al ) and 70 times lower than its peak flaring flux (jAharonian et al.l 120071 ). The positive 
excess variance a^ s , in dicating variability, a llows a fractional root mean square (rms) of 



F var V HE = 23 ± 3% (see lVaughan et al.ll2003l for definitions of a^s an d F VF , r ) to be derived. 



which is 3 times less than the high state flaring variability reported by lAharonian et al. 



(120071 ) . A spectrum was obtained for each night when possible, otherwise two or three nights 
were combined. No indication of spectral variability was found during those observations, 
with a limit on the nightly index variations of Ar < 0.2. 

The Fermi light curve shows the photon fluxes for the high energy (HE) range, 0.2- 
300 GeV, and the photon spectral indices for each interval. Each bin is the result of a power 
law fit, using the background values found on the overall time-averaged fit, and centered 
on the H.E.S.S. observations. The light curve fit to a constant has a x 2 probability of 
p{x 2 ) = 0.95, clearly consistent with a constant flux. The normalized excess variance of 
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—0.16 ± 0.09 sets a 90% c onfidence level limit of -F V ar,HE < 20% on the fractional variance 
jFeldman fc Cousins! ^998^ . 



The X-ray light curve, derived from spectral fits of the nightly RXTE (and Swift) data 
sets, shows flux doubling episodes on time scales of days, similar to the optical and VHE 
measurements. The lowest fluxes of ~ 3-6 x 10~ 1 1 ergcm~ 2 s~ 1 are at the same level as 
those seen in the low state (jAharonian et al.l l2005bl ) but with larger fluctuations, F var) x = 



35 ± 0.05%. The time history of the fitted spectral indices in Fig. [T] show clearly that the 
X-ray spectrum hardens significantly, AT X 0.5, as the 2-10 keV flux increases. 



The ATOM fluxes are ~ 5 times higher than the low state found in lAharonian et al. 



(j2005bl ). but the V^-band magnitudes reported here are in the ra nge 12.7-13 which is w ell on 
the lower side of the measurements of PKS 2155—304 reported by lFoschini et al.l (120081 ) when 
the source was quoted to be in a low state with l/-band ma gnitudes in the range 1 2-12.7. 
The host galaxy flux is estimated to be ~ 10 -11 ergs cm -2 s _1 (jKotilainen et al.lll998l ). hence 
most of the optical flux can be attributed to the central AGN. The average fractional rms 
over all bands is -F var opt ~ 8 ± 0.5%. The B — R lightcurve is compatible with a constant, 
p(x 2 ) = 0.66, indicating little or no optical spectral variability. 



3. Discussion 

The two-component broad band spectra of high energy-peaked BL L ac objects (HBLs) 



are ty pically modeled with synchrotron self-Compton (SSC) scenarios (e.g-. lBand fc Grindlay 



19851 ). Despite the simplicity of these models, they have been successful in reproducing 
many blazar SEDs and make definite predictions for the flux and spectral variability that 
should be seen in the two components. In particular, for typical parameters, the electrons 
responsible for the X-ray emission also produce the VHE emission; and if the underlying 
particle distributions were to vary, the resulting flux and spectral changes in the VHE band 
should be related to variations in the X-rays. In fact, for the July 2006 flare, a non-linear 
relationship was seen between the X-ray and VHE band s, though the observed variability 



patterns do not quite fit the simple SSC model in detail (ICostamante et al.l 120081 ) . 



In Fig. [21 we overlay a model SSC calculation that roughly fits the time-averaged SED. 
The electron distribution model parameters, a three-component power-law with indices p^ = 
1.3, pi = 3.2, p 2 = 4.3 (dn/d'j oc / y~ Pi ), minimal and maximal Lorentz factors 7 m j n = 1 
and jmax = 10 6 ' 5 , break electron Lorentz factors 71 = 1.4 x 10 4 , 72 = 2.3 x 10 5 , and total 
electron number N tot = 6.8 x 10 51 , have been set to reproduce the shape of the lower energy 
component of the SED. The overall SED is then adjusted with the remaining parameters: 



- 11 - 



radius of the emitting region in the comoving frame, R = 1.5 x 10 17 cm; bulk Doppler factor, 
5 = 32; magnetic field, B = 0.018 G. Even though we regard this fit as a "straw-man" 
model, it is perhaps reassuring that the joi nt Fermi-H.E.S.S. time-ave raged spectra can be 



reasonably well-described as SSC emission. iKatarzynski. et al.l (120081 ) found similar values 



for R, B and 6 in their SSC description of a steady large jet component in the SED of 
PKS 2155-304. 

Some features of this model calculation are particularly noteworthy. The electrons that 
produce the synchrotron X-ray emission have Lorentz factors > 72. When the power-law 
component for those electrons is omitted from the calculation, the dot-dashed curve in Fig. [2] 
results. For this particular set of parameters, the electrons that produce the X-rays have 
higher energies than the electrons that produce the VHE emission. Furthermore, the lack of 
a significant impact on the shape of the SSC component when those electrons are removed 
indicates that Klein-Nishina effects suppress any significant contribution by those electrons 
to the emission at ~ TeV energies. 

These features of this calculation allow that there need not be a correlation between 
the X-ray and VHE fluxes; and in fact, this is what is observed. In contrast with the July 
2006 flare, we do not find any evidence of flux correlation between the X-ray and H.E.S.S. 
bands with a Pearson's r of 0.12 ±0.1 between these bands. Furthermore, the 2-10 keV X- 
ray spectra show spectral variability consistent with an underlying electron distribution for 
which the cooling time scales are of order the flux variability time scales, i.e., the spectra are 
softer when the flux is lower, with changes in photon index of Ar^. « 0.5 (Fig. [T]); whereas 
the VHE emission shows no evidence for significant spectral variability despite flux variations 
of a factor of 2. Since radiative cooling time scales vary inversely with electron energy, this 
supports the conclusion that the electrons responsible for the synchrotron emission in the X- 
ray band have higher energies than the electrons that produce the inverse- Compton emission 
in the VHE range, assuming they are part of the same overall non-thermal distribution. 

Even though this all fits in with our straw-man SED calculation, the variability patterns 
in the optical, X-ray, HE and VHE bands suggest a much more complex situation. In 
the absence of spectral variability, the mechanisms that would produce the observed flux 
variability in the VHE band are rather constrained. Increases in flux could be driven by 
injection of particles with a constant spectral shape, and decreases in flux could be caused 
by particle escape from the emitting region or by expansion ( "adiabatic" ) losses, assuming 
those latter two processes can operate independent of particle energy. 

However, since the electrons that produce the VHE emission must be in the weak radia- 
tive cooling regime, a more natural mechanism for the flux variability would be that changes 
in the seed photon density are driving the variability. Comparing the daily flux values in 
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the optical and the VHE bands, we find indications of fairly strong correlations that sug- 
gest that the optical emission provides the target photons for the IC emission. In the B, 
V, and R bands, the correlations with the H.E.S.S. fluxes have Pearson's r values in the 
range 0.77-0.86 with uncertainties < 0.09. This correlated behavior is readily apparent in 
the light curves shown in Fig. [TJ and these results provide the first quantitative evidence of 
correlated variability between the optical and VHE bands on these time scales for an HBL@ 
Confirmation of this behavior, not only from this source but also from other VHE emitting 
blazars in a low state, would provide important constraints on emission models for these 
objects. 

In the context of a single-zone SSC model, we would expect that any flux variability 
in the optical bands should also appear as variability in the Fermi-LAT energy range. To 
illustrate this, we plot, as the dashed curve in Fig. [2] the SED that results if we omit contri- 
butions from electrons with energies > 71. For the original model parameters, the electrons 
that produce the optical-soft X-ray emission also produce the bulk of the IC component, 
including the HE and VHE emission. Since we do not find any indication of a correlation 
between the optical and HE fluxes, this suggests that the optical emission may arise from a 
separate population of electrons than those responsible for the HE and VHE emission. If so, 
then these electrons probably also occupy a distinct physical region with different physical 
parameters (magnetic field, size scale, bulk Lorentz factor). Multizone SSC models of this 
kind have already been proposed to accoun t for the "orphan" 7-ray flare in 1ES 1959+650 



during May 2002 flKrawczvnski et al. 112004 ). 



Although the 0.2-300 GeV photon fluxes measured by Fermi are consistent with being 
constant, we find more significant variations of the photon spectral index in the daily analyses 
{p{x 2 ) = 0.19). The fitted values range from fairly soft, T = 2.7 ± 0.7, to extremely hard, 
r = 1.1 ±0.4. 

These values, along with the constant, intrinsic VHE index of Tvhe ~ 2.5 derived from 
the H.E.S.S. data, imply spectral breaks between the HE and VHE bands of Ar as large as 
1.4. Very sharp spectral breaks (Ar > 1) would require rather narrow electron distributions 
and would therefore pose difficulties in fitting a broad lower energy component in the context 
of a single-zone model. Interestingly, we find a significant anticorrelation between the nightly 
X-ray fluxes and the Fermi-LAT spectral indices of rxr = —0.80 ± 0.15. A fit to a linear 
model is preferred over a constant at the 2.6<r level, with a slope of —0.14 ± 0.05. If the 
electrons that produce the X-rays are at higher energies than those that produce the TeV 



Donnarumma et al.l (|2008f ) mention possible correlated variability in the recent June 2008 flare of Mrk 421 



in a high state. 
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emission, the cause for such a correlation would be difficult to understand. An important 
caveat in considering these results is that the Fermi coverage for PKS 2155—304 was relatively 
uniform over each 24 hour period, whereas the optical, X-ray, and VHE observations were 
restricted to 4-6 hour intervals each night. Hence, the Fermi observations are not strictly 
simultaneous with the other measurements, so it is possible that some of the observed HE 
spectral variability occurred outside of the nightly observing windows. 

As the first multiwavelength campaign of an HBL that includes Fermi and an ACT 
instrument, these observations have yielded results that strongly challenge the standard 
models for these sources. Having caught PKS 2155—304 in a low state, we see that its spectral 
and variability properties are significantly different than its flaring, high state behavior. The 
variability patterns, in particular, defy easy explanation by the usual SSC models and should 
provide valuable constraints for models that attempt to describe the emission mechanisms 
in blazar jets. 
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Fig. 1. — Light curves from (top to bottom): H.E.S.S., Fermi, RXTE/Swift, and ATOM. 
The Fermi and RXTE/Swift panels also show the spectral index measurements (red) for each 
night. Vertical bars show statistical errors only. Horizontal bars represent the integration 
time and are apparent only for the RXTE and Fermi data. The ATOM bands are B (blue 
circles), V (green squares) and R (red squares). 
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Fig. 2. — The SED of PKS 2155—304. The red butterfly is the Fermi spectrum restricted 
to the MJD 54704-54715 period, while the black butterfly covers MJD 54682-54743. As a 
cross check of the fit robustness, the differential flux was estimated in 8 limited energy bins 
by a power law fit (black circles) and are found to be consistent within la of the global fit, 
including a clear spectral break at ~ 1 GeV. The gray points are archival NED data, and 
the two gray butterflies are EGRET measurements. The solid line is a 1-zone SSC model. 
The dashed and the dot-dashed lines are the same model without electrons above 71 and 
72, respectiv ely. The VHE pa r t is ab sorbed with the P0.45 extragalactic background model 
described in lAharonian et all (j2006al ) . 



